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ABSTRACT 

This paper is the 4th in a series describing the latest additions to the BaSTI stellar evolution 
database, which consists of a large set of homogeneous models and tools for population synthesis 
studies, covering ages between 30 Myr and ~20 Gyr and 11 values of Z (total metallicity). Here we 
present a new set of low and high resolution synthetic spectra based on the BaSTI stellar models, 
covering a large range of simple stellar populations (SSPs) for both scaled solar and a-enhanced 
metal mixtures. This enables a completely consistent study of the photometric and spectroscopic 
properties of both resolved and unresolved stellar populations, and allows us to make detailed 
tests on specific factors which can affect their integrated properties. Our low resolution spectra 
are suitable for deriving broadband magnitudes and colors in any photometric system. These 
spectra cover the full wavelength range (9-160000nm) and include all evolutionary stages up to 
the end of AGB evolution. Our high resolution spectra are suitable for studying the behaviour 
of line indices and we have tested them against a large sample of Galactic globular clusters. We 
find that the range of ages, iron abundances [Fe/H], and degree of a-enhancement predicted by 
the models matches the observed values very well. We have also tested the global consistency 
of the BaSTI models by making detailed comparisons between ages and metallicities derived 
from isochrone fitting to observed colour-magnitude diagrams, and from line index strengths, for 
the Galactic globular cluster 47 Tuc and the open cluster M67. For 47 Tuc we find reasonable 
agreement between the 2 methods, within the estimated errors. From the comparison with M67 
we find non-negligible effects on derived line indices caused by statistical fluctuations, which are 
a result of the specific method used to populate an isochrone and assign appropriate spectra to 
individual stars. 

Subject headings: stars: evolution — galaxies: evolution — galaxies: stellar content 

1. Introduction both test the theory of stellar evolution, and in- 

vestigate the formation and evolution of galaxies, 
Large grids of stellar models and isochrones are Qne of the major open problems of mo dern astro- 

necessary tools to interpret photometric and spec- physics 

troscopic observations of resolved and unresolved _ ,, . , . 

„ , . . „ lo this purpose, we started m 2004 a major 

stellar populations, this, m turn, allows one to . , , , , 

project aimed at creating a large and nomoge- 
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neous database of stellar evolution models and 
isochrones (BaSTI - a Bag of Stellar Tracks and 
Isochrones) that cover a large chemical composi- 
tion range relevant to stellar populations in galax- 
ies of various morphological types. In addition, 
BaSTI includes the option of choosing between 
different treatments of core convection and stel- 
lar mass loss employed in the model calculations. 
In its present form, BaSTI contains a grid of mod- 
els that cover stellar populations with ages rang- 
ing between 30 Myr and ^20 Gyr, which include 
all evolutionary phases from the Main Sequence 
(MS) to the end of the Asymptotic Giant Branch 
(AGB) evolution or carbon ignition, depending 
on the value of stellar mass. They are calcu- 
lated for two metal mixtures (scaled solar and a- 
enhanced) and 11 values of Z (total metallicity), 
each with tw o choices for t he Reimers mass loss 
parameter r\ ( Reimers 19751 ). and two mixing pre- 
scriptions during the Main Sequence (without and 
with overshooting from the Schwarzschild bound- 
ary). Broadband magnitudes and colors for the 
BaSTI models and isochrones are provided in vari- 
ous photometric systems (Johnson-Cousins, Sloan, 
Stromgren, Walraven, ACS /H ST), making use of 
bolometric corrections derived from an updated 
set of model atmosphere calculations, for element 
mixtures consistent with those employed in the 
stellar evolution computation^- 

In addition to stellar models and isochrones, 
BaSTI also provides a series of Web tools that 
enable an interactive access to the database and 
make it possible to compute user-specified evolu- 
tionary tracks, isochrones, stellar luminosity func- 
tions (star counts in a stellar population as a func- 
tion of magnitude) plus synthetic color-magnitude 
diagrams (CMDs) for arbitrary star formation his- 
tories (SFHs). The results of this maj o r effort have 
been published in IPietrinferni et al.l (120041. 12 006. 
herein after Papers I and II) and [Cordier et al.1 
(|2007l hereinafter Paper III). All Web tools, mod- 
els and isochrones can be found at the BaSTI offi- 
cial website |http -J/193. 204 . 1 . 62/index . html 



amp les, the tests presen t ed in Papers I, II an d 
III; ISalaris et all l2007t iTomasella et al.l [2008). 



The database has been extensively tested 
against observations of local stellar populations 
and eclipsing binary systems (see, for some ex- 



X A11 the results presented in this work are based on the 2008 
release of the BaSTI archive. We refer to the official BaSTI 
website for more details on this new release. 



In its present form BaSTI can be used to in- 
vestigate resolved stellar populations, and in- 
deed it has been widely employed in studies of 
Galactic and extragalactic resol v ed star clusters 



see, 



e.g. De Angeli et al. 20051 : Villanova et al 



20071: iMackev & Brobv Nielsen] l2007l for just 



a few examples) and in the determination of 
SFH and chemical enrichment histo ry of resolved 
Loca l Group galaxies (see 



2004; 



Barker et al 



■R- 



Gallart et al 



20071: iGullieuszik et al.l 120071: 



Carrera et al.ll2008h . 

Integrated colors and magnitudes for sim- 
ple (single-age, single metallicity) stellar popu- 
lations (SSPs) can be easily determined from the 
isochrones. The only additional pieces of informa- 
tion needed are an Initial Mass Function (IMF) 
and the relationship between the initial and fi- 
nal (remnant) stellar masses. Indeed, some in- 
tegrated colors from BaSTI models have already 
been tested against empi rical constra i nts in Pa- 
pers II and III. M oreov er, I James et al. I (|2006l) and 
ISalaris fc CassisJ (|2007l ) have already employed in- 
tegrated colors from BaSTI isochrones to address 
issues related to the SFH of unresolved dwarf and 
elliptical galaxies, and extragalactic globular clus- 
ter ages. 

A further development of BaSTI, which we 
present in this paper, is the inclusion of integrated 
spectra (including spectral line indices) for SSPs 
spanning the same parameter space - in terms of 
ages, metallicities and heavy element mixtures - 
covered by the isochrones of Papers I and II. Cor- 
responding quantities for composite stellar popu- 
lations with an arbitrary SFH can then be easily 
determined by integration over a range of discrete 
SSPs. 

We present both low and high resolution ver- 
sions suited to different purposes. Our low resolu- 
tion spectra cover all evolutionary stages and are 
based on libraries of (predominantly) theoretical 
stellar spectra, consistent with those used to pro- 
duce the bolometric corrections employed in the 
isochrone calculations. This ensures that our SSP 
integrated broadband magnitudes and colors ob- 
tained from either adding up the broadband fluxes 
of the individual stars populating the appropri- 
ate isochrone, or from their composite integrated 
spectrum, will be exactly the same (see test in sec- 
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tion !2.ip . Our high resolution spectra are suitable 
for deriving line strengths and making detailed 
tests of their behaviour as a function of age and 
chemical composition, and also horizontal branch 
morphology and sampling in the case of resolved 
stellar clusters (see tests in sections 13.11 and I3.2[) . 

With the inclusion of predictions of both photo- 
metric and spectroscopic integrated properties of 
SSPs, BaSTI will provide a well tested, fully ho- 
mogeneous and up-to-date set of theoretical tools 
that enable us to investigate self-consistently both 
resolved and unresolved stellar populations. We 
present a global test of this self-consistency and 
reliability of the database, by considering both 
the observed CMD and integrated spectrum of the 
well-studied clusters 47 Tuc and M67, to inter- 
compare ages and metallicities obtained from our 
high-resolution spectra with the results from fit- 
ting theoretical isochrones to the observed CMD, 
and with independent chemical composition esti- 
mates from spectroscopy of individual stars. 

The structure of the paper is as follows. A 
description of the method used to create our theo- 
retical integrated spectra is given in Section [2j 
including separate detailed descriptions of our 
low and high resolution spectra in Sections 12.11 
and 12.21 respectively. This is followed in Sec- 
tion [3] by a comparison of our derived line in- 
dices with the Galact ic globular cluster sample of 
Schiavon et al. ( 2005f) . A global consistency test 
of our isochrones, integrated colors and spectra us- 
ing photometric and spectroscopic data for 47 Tuc 
and M67 is presented in Sections 13.11 and 13.21 A 
brief summary in Section [4] closes the paper. 

2. Creating integrated spectra: Method 

Before an integrated spectrum can be created, 
the SSP (or composite population) itself must be 
created by 'populating' the relevant isochrone (or 
combination of isochrones). BaSTI isochrones are 
defined in terms of total metallicity Z (with cor- 
responding [Fe/H] depending on whether scaled 
solar or a-enhanced models are being used) and 
age. Each isochrone consists of 2250 discrete evo- 
lutionary points, each of which is defined in terms 
of effective temperature T e ff, luminosity L, and 
mass M (from which logg can also be derived). 
An isochrone i s populated by applying the IMF by 
Krouoal (|200lh with an appropriately chosen nor- 



malization constant - the value of this constant is 
given by M t = J^ u -0(M)MdM where M t is the to- 
tal mass of stars formed in the initial burst of star 
formation that originated the SSP. The upper and 
lower integration limits correspond to the lowest 
and highest mass stars formed during the burst, 
and are fixed to 0.1 and 100M Q , respectively. We 
choose to set M t equal to IMq. We note here 
that the lower mass limit of the BaSTI database 
is 0.5 M & , however objects with masses below this 
limit contribute negligibly to th e integrated mag- 
nitude s and colors, as verified in Salaris k, Cassisil 
(l2007h by implementing the very low mass star 
models bv lCassisi et al. (2001), that extend down 
to ~0.1Mq. Hence, stars of < 0.5 M Q are in- 
cluded in the total mass and are accounted for in 
the normalization constant, but their photomet- 
ric properties are not included in the integrated 
spectrum. 

Once an SSP has been created, each evolution- 
ary point in a particular simulation is assigned a 
spectrum by interpolating linearly in [Fe/H], T e ff 
and log<? amongst spectra in the appropriate spec- 
tral library (see Sections 12.11 and 12.21 for choices 
of spectral libraries). Each spectrum is scaled by 
the stellar surface area, via the stellar radius de- 
rived from the evolutionary point on the isochrone, 
and then weighted appropriately, according to the 
number of stars predicted by the normalized IMF, 
as described above. The 2250 individual spectra 
are then summed to obtain the final integrated 
spectrum. 

2.1. Integrated spectra: low resolution 

Our low resolution integrated spectra incorpo- 
rate all evolutionary stages covered by the BaSTI 
isochrones, as detailed in Section [TJ This yields 
a full spectral energy distribution (SED) suitable 
for calculating broad-band colors and, for exam- 
ple, deriving K-corrections. 

The majority of the low resolution spectra 
used i n this study are from the Castelli fc Kuruczl 
(|2003h data set@, based on the ATLAS9 model at- 
mospheres. These spectra cover effective temper- 
atures, Teff, from 3500K up to 50000K and logg 
from 0.0 dex to 5.0 dcx. Metallicities range from 
[Fe/H]= —2.5 to +0.5 for both scaled solar and 
a-enhanced mixtures - the level of a-enhanccmcnt 



2 available at http://kurucz.harvard.edu/grids.html 
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being fixed at [a/Fe]=0.4, the same as that used in 
the BaSTI stellar models. These spectra cover the 
full wavelength range (9-160000nm) and sampling 
varies from lnm at 300nm to lOnm at 3000nm. 

For stars with T e tf < 3500K (except for car- 
bon stars in the AGB phase - see below) we sup- 
plement the Castelli/Kurucz spectra with those 
from the BaSeL 3.1 (WLBC 99) spectral library 
(jWestera et al.ll2002ft . The WLBC99 spectra have 
the same wavelength coverage and sampling as 
the Castelli/Kurucz ones and are provided for 
scaled solar abundances only, in the range —2.0 < 
[Fc/H] < +0.5. Since we are only using the 
WLBC99 spectra for the coolest stars (mostly 
RGB, and some AGB), the lack of a-enhanced 
spectra should not be a problem as these stars 
only contribute significant flux in the reddest part 
of the spectrum where it is known that the colors 
are gene rally insensitive to the metal distribu tion 
(see e.g. lAlonso et ail 1 19991: ICassisi et "aLll2004h . 

For the AGB carbon stars we use the aver- 



aged AGB spectra of lLancon &; Mouhcind (|2002l ). 
These (empirical) spectra cover 510-2490nm, 
with linear extrapolation to zero flux at 350 and 
5000nm, at a resolution of 0.5nm (5A), which 
we rebin to match the Castelli/Kurucz sampling. 
They are defined in terms of T e ff only, as no 
metallicity information is available. 

Some representative examples of our low reso- 
lution integrated spectra are shown in Figure [1] 
Since the low resolution spectra have full wave- 
length coverage and cover all evolutionary points 
on the isochrones, they are useful for identifying 
the flux contribution from specific stars or evolu- 
tionary stages. This is important when we assess 
the validity of line strengths derived from the high 
resolution spectra, for which there is a low tem- 
perature cut-off in the spectral library used (see 
section l2~2l) . 



2.1.1. Integrated colors - behaviour and compar- 
isons 

As a sanity check, we verified the consistency 
of the colors derived from our low resolution in- 
tegrated spectra with those predicted from the 
isochrones by the BaSTI population synthesis Web 
tool. Colors were derived from the spectra by con- 
volving the integrated spectra with broadband fil- 
ter profiles and normalising to a synthetic Vega 



spectrum ( Castelli fc Kuruc3 [19941 . Colors are 
obtained from an isochrone by adding up the flux 
contributions from each evolutionary point in the 
relevant broadband filter. For a given age t, metal- 
licity Z, plus an IMF ip(M), we simply integrate 
the flux in the generic broadband filter A along the 
representative isochrone, i.e. 



M\ = -2.5 log 




^(M)10-°- 4M ^ M )dM 



We performed this test for a typical globular clus- 
ter (Z=0.001, i=10Gyr) and a typical open cluster 
(Z=0.0198, t=3.5Gyr), from (U - B) through to 
(J — K). For the solar metallicity cluster, all col- 
ors are consistent within < 0.015 mag whilst for 
the lower metallicity simulation the agreement is 
even better (within < O.Olmag). 

BaSTI isochrones currently provide magnitudes 
for 'standard' UBVRIJHKL filters. Since we 
have shown that the resulting colors are consis- 
tent with those derived from our integrated spec- 
tra, magnitudes and colors in other photometric 
systems can reliably be derived by convolving the 
spectra with the required filter response curves, 
and using the appropriate normalizing spectrum 
(as described above). 

As an example of our integrated SSP color pre- 
dictions, we display in Figure[2]the run of (U — B), 
(B - V) and (J - K) as a function of the SSP 
age, for the scaled solar and a-enhanced mixtures. 
Scaled solar models are shown for [Fe/H]=+0.06, 
—0.66 and —1.27, and are matched to the a- 
enhanced models with almost exactly the same 
[Fe/H] values ([Fe/H]=+0.05, -0.70 and -1.31). 
Differences between the results for the two mix- 
tures are overall reasonably small, especially for 
ages above 1 Gyr. For these old populations 
the largest color variations at a fixed age are 
^0.05 mag for (J — K) and (V — I), in the sense 
of the a-enhanced SSPs (that have a higher Z at 
fixed [Fe/H]) being redder. The (U — B) colors 
at fixed [Fe/H] for the two metal distributions are 
overall more similar than the case of (J — K) and 
(V-I). 

We can also compare our integrated colors 
to the predictions from the very recent models 
by ICoelho et all (|2007h . They computed inte- 
grated spectra and colors for a range of SSPs 
with scaled solar and a-enhanccd metal mixtures 
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using a method which is qualitatively similar to 
ours. Their choice for the a-enhancement is essen- 
tially equal to ours, but the stellar isochrones and 
spectra employed in their modelling arc different, 
and they also cover only a restricted metallic- 
ity range at present, with [Fe/H] between —0.5 
and +0.2. Their [Fe/H] gridpoint closest to 
ours is [Fe/H]=0.0 for both scaled-solar and a- 
enhanced mixtures, and so we can compare these 
to our scaled solar [Fe/H]=+0.06 and a-enhanced 
[Fe/H] =+0.05 models. Figure [3] displays the 
(U-B), (B-V) and (V-I) integrated colors for 
the common age range. Our colors are systemati- 
cally redder for both a-enhanced and scaled-solar 
mixtures, except for the (V — I) colors, where the 
reverse is true. Typical differences are of the order 
of ~0.1 mag. In section [2. 2. 11 which compares the 
behaviour of line ind ices predicted by t he 2 sets 



of models (ours and ICoelho et al.ll2007T) . we dis- 



cuss some factors which may contribute to these 
differences. 

2.2. Integrated spectra: high resolution 

Our high resolution integrated spectra are suit- 
able for deriving line strengths (e.g. of Lick-style 
indices) and studying their behaviour as a func- 
tion of age, metallicity and a-enhancement (see 
section l2~.2.ip . We also use them to make some pre- 
liminary tests of the effects of, for example, hor- 
izontal branch morphology, and statistical sam- 
pling in resolved clusters (see section [3]). They 
are const ructed using the libra ry of s ynthetic spec- 
tra fro m iMunari et all (|2005l ) . The IMunari et all 
(|2005l ) spectr a were compu ted using the SYN- 
THE code byE uruca (|1993l ) and the input model 
atmospheres are the same as those employed to 
calculate the l ow-resolution spectral li brary de- 
scribed above ( Castelli k. Kuruc3 2003h - hence 
they have similar metallicity, T e // and logg cov- 
erage for bo th scaled sola r and a-enhanced ver- 
sionaj. The Munari et al. (|2005l) high resolution 
spectra take into account the effects of several 
molecules, including C 2 , CN, CO, CH, NH, SiH, 
SiO, MgH, OH, TiO and H 2 0. The source of 
atomic and molecular data is the Kurucz a tomic 
and molecular line-list (see, e.g. Kuruczfl992l ) with 
the exception of TiO and H 2 data, for which 



the line lists of, respect i vely, 



Partridge fc Schwenkd (|1997I ) were adopted. The 



Schwenkd (|l998f) and 



wavelength range of the Munari spectra is 2500- 
10500A and the authors provide several different 
resolutions - we used those with a uniform disper- 
sion of lA/pix here. 

The referee suggested that we s hould investi- 



gate t he comparison between the IMunari et al 



(|2005l ) theoretical stellar spectra and their empir- 
ical counterparts, particularly for cool stars, for 
which Arcturus is the best studied example. In 
fact, our adopted high-resolution spectral library 
has bee n extensively tested again st observed spec- 
tra by iMartins fc Coelhol (|2007h . In particular, 
they measured 35 spectral indices defined in the 
literature on three synthetic spe ctral databases 
(including the IMunari et al. 20051 library adopted 
in this work) and compared them with the cor- 
responding values measured on th ree empirical 



spect ral librari es, namely Indo-US (IValdes et al . 
2004) , MI LES IsSichez-Blazquez et al.l l2006l) and 



ELODIE (jPrugniel fc Soubiranll200l[ ). The reader 
is referred to this very informative paper for de- 
tails about the comparisons and results. As a gen- 
eral conclusion, all three sets of synthetic spectra 
tend to show the largest discrepancies with the em- 
pirical counterparts at the coolest temperatures. 
However, the quantitative and qualitative trends 
of these differences depend on the specific empiri- 
cal dataset used for comparison. 

We now consider briefly the H/3 and Fe5406 
Lick-style indices, which we will be using exten- 
sively in the rest of this paper (see Figs. 29, 35 and 
Tables CI, C2 and C3 in Martins & Coelho 2007). 
For H/3 the best agr eement, as given by the a dev 
parameter defined in IMartins fc Coelhol (|2007ft . is 
with the Indo-Us library for T eff >7000 K (high 
temperatures) and T e // <4500 K (low tempera- 
tures), whereas the best agreement at intermedi- 
ate temperatures is with the MILES library (the 
reader is warned that, judging from Tables CI, 
C2 and C3, some of the adev values in their Figs. 
29 and 35 have been mistyped). Fo r a solar-like 



dwarf star, the IMunari et alJ (|2005l ) spectra re- 



3 see http://archives.pd.astro.it/2500-10500/ for details of 



coverage 



produce well the observed value, whereas for an 
Arcturus-like giant the theoretical H/3 index is 
lower than observations. In general, for all three 
temperature ranges there is no clear systematic 
offset between theory and observations, but a non- 
negligible spread in the difference between them. 
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In the case of Fe5406 the best agreement is 
with the MILES database at intermediate and 
high temperatures, with ELODIE being only 
marginally more in agreement than MILES at 
low temperatures. The Indo-US library displays 
a clear systematic difference at low temperatures 
(theoretical indices being larger) and, to a lesser 
extent, also in the intermediate regime. These 
offsets are much less evident in comparisons with 
the ELODIE database. For a solar-like dwarf star 
our adopted spectra tend to slightly overpredict 
the observed value, whereas for an Arcturus-hke 
giant the theoretical index is in agreement with 
observations from the ELODIE library, slightly 
overpredicted compared to the Indo-US database, 
and massively overpredicted compared to MILES. 

The reasons for the discrepancies with observed 
spectra are certainly at least partly due to inad- 
equaci es in the theoretica l spect ra, but, as noted 
also bv lMartins &: Coelhol ( 2007 ). errors in the at- 
mospheric parameters of the observed stars may 
also play a crucial part, especially in explaining 
the different trends with observations coming from 
the different empirical libraries. 

To this purpose, we made our own test on 
Arcturus, using an empirical spectrum from the 
ELODIE spectral library. We used their 0.2A res- 
olution spectrum, which we degrad ed to lA reso- 



lution to match the resolution of the Munari et al 



(|2005l ) library. We compared this empirical spec- 
trum with a theoretical one which we created by 
interpolating amongst t he Munari spectra, us ing 
stellar parameters from IPeterson et aD <|l993h - 
namely [Fe/H]=-0.5 (±0.1), T e// =4300 (±30 K) 
and logg=1.5 (±0.15). W e used an q-enhanced 
theoretical spectrum since IPeterson et al. (1993) 
found a-elements enhanced by ~0.3 dex, although 
it should be noted that the theoretical spectra are 
fixed at [a/Fe]=0.4. 

Measuring H/3 and Fe5406 line strengths di- 
rectly on these spectra (see section 12.2.11 for more 
details on this), we found that H/3 is under- 
predicted in the synthetic spectrum by 0.25A 
(measuring 0.6lA, whereas the empirical spec- 
trum gives 0.86A), whilst the Fe5406 line is slightly 
over-predicted (2.27A cf. 2.18A). We noted that 
the Mgb line is also overpredicted by a simi- 
lar fraction. An underpredicted value for the 
H/3 index is most likely explained by the lack 
of non-LTE and/or chromospheric contributions 



and /or inadequate treatment of convection ( see , 
e.g. iMartins fc Coelhol l2007t iKorn et al.ll2005h in 
the theoretical spectra. However, the effect of 
uncertainties in the atmospheric parameters de- 
termined from empirical spectra may play also a 
crucial part in explaining these differences. As a 
test, we modified T e ff and logg in the synthetic 
Arcturus comparison spectrum to see the effect on 
derived indices. We found that by increasing T e f / 
by 70K and decreasing logg by 0.25 dex (~ 2a 
change in both parameters) the H/3 line strength 
increased to 0.83A and Fe5406 decreased to 2.12A, 
both within 3% of the empirical values. The Mg6 
line strength also decreased to come into closer 
agreement with the empirical spectrum. 

The full integrated SED derived from the high 
resolution spectra differs fro m that of the l o w res - 
olution spectra because the Munari et al. ( 2005h 



spectra do not include temperatures cooler than 
T e ff = 3500K (the low temperature limit of the 
ATLAS9 model atmospheres) and hence the full 
SED differs from that of the low resolution spec- 
tra, particularly longward of- 6000A. This means 
that our high resolution integrated spectra are, in 
general, not suitable for deriving broad band col- 
ors and hence we use the low resolution versions 
for this. Figure [1] shows some representative ex- 
amples of our high resolution integrated spectra. 

We decided not to attempt to include (or 
create) lower temperature high resolution spec- 
tra because we are concerned about possible in- 
consistencies between stellar atmosphere mod- 
els from different sources. Figure [4] shows a 
co mparison betwe e n lin e strengths derived from 



the iMunari et al.1 (|2005j ) spectra (based on AT- 
LAS9 model atmospheres) and corresponding 
spectra from the MARCS databas^l - we com- 
pare strengths of the Hf3, Mgb and Fe5406 lines 
since we will be using these extensively in our 
more detailed analyses later in this paper. It can 
be seen that for all 3 lines there are significant off- 
sets, and that the behaviour becomes increasingly 
discrepant between the 2 models as the tempera- 
ture decreases. We note that these comparisons 
were done using the MARCS models available as 
of May 2008. New low temperat ure MARCS mod- 



els ar e currently being produced ( G ustafsson et al 



2008) - we will investigate these as they become 



4 http:/ /marcs. astro. uu.se/ 
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available and may include them in our integrated 
spectra at a later date. We note here that this low 
temperature regime also includes cool AGB stars, 
for which no realistic models currently exist. 

The effect of not including stars with T e tf < 
3500K in the high resolution integrated spectra is 
illustrated in Figure [5] Here we use the low res- 
olution spectra to separate out the flux contribu- 
tion from these stars, for two representative SSPs, 
and plot the percentage of the total flux coming 
from T e ff < 3500K as a function of wavelength. 
It can clearly be seen that the flux contribution 
from these stars only becomes significant at wave- 
lengths longer than the bandpasses of the main 
Lick-style indices, e.g. H/3, Mg& and various Fe 
lines f see 12.2. ip . 

It is expected that high metallicity SSPs will 
be most affected by this low temperature cut-off 
in the spectral library, since a larger portion of 
their isochrones fall in the T e ff < 3500K regime 
than for low metallicity SSPs. As a test of the 
maximum likely effect on various diagnostic lines 
(see next section for definitions), we mimicked 
the effect of a T eff = 3500K cut on a Z=0.04 
(i.e. ~ twice solar), 10 Gyr SSP by taking a 
10 Gyr, Z=0.001 (~ ^ solar) isochrone and im- 
posing a temperature cut at T e ff =4I30K. This 
temperature cut corresponds to approximately the 
same evolutionary point on the isochrone as a 
T eff = 3500K cut for the higher metallicity SSP. 
We then created an integrated spectrum for this 
'truncated' Z=0.00I SSP and compared the de- 
rived line indices with those derived from the in- 
tegrated spectrum for the full SSP at the same 
age and metallicity. For the 'truncated' spectrum, 
the H/3 index (the main age indicator - see next 
section) was ~ 0.1 higher than that of the full 
spectrum. This corresponds to an age difference 
of ~ 1 Gyr, and goes in the sense that the trun- 
cated spectrum looks too young. There are vari- 
ous metallicity indicators (see next section) - for 
the main [Fe/H] indicator we use here, Fe5406, the 
line strength is slightly lower than for the full spec- 
trum, which would result in an underestimate of 
the true [Fe/H] by ~ 0.05 dex. We repeated this 
test for a 3 Gyr model and found similar absolute 
offsets, which would result in slightly smaller ef- 
fects on the inferred age, at the level of ~ 0.5 Gyr. 

It should be noted that the issue of 'missing 
flux' at low temperatures, and its associated ef- 



fect on line indices, is not a problem specifically 
associated with the use of theoretical spectral li- 
braries. It is also very relevant to any population 
synthesis techniques which rely on empirical spec- 
tra, since there are very few local stars at these 
low temperatures, particularly at low metallicity. 

2.2.1. Line indices: behaviour as a function of 
age and metallicity 

All the line indices discussed in this pa- 
per are defined by the bandpasses tabulated in 



iTrager et al.1 (|1998l i.e. the 21 classical Lick/IDS 
indices). The quoted line strengths, and hence 
those displayed in plots, were obtained using the 
LECTOR programme by A. Vazdeki^l and are 
those directly measured on the spectra, i.e. they 
are not transformed onto the Lick system. 

Before discussing the behaviour of various line 
indices it is worth remembering the distinction be- 
tween iron abundance, quantified here as [Fe/H], 
and total metallicity, Z . For scaled solar models, 
all metal abundances scale as for the sun, whilst 
in a-enhanced models the a element abundances 
are increased relative to Fe. Hence at fixed Z, the 
a-enhanced models have a lower [Fe/H] than the 
corresponding scaled solar models. As an example, 
the BaSTI scaled solar models with Z=0.0I98 (the 
solar value) have [Fe/H] =+0.06, whilst the corre- 
sponding a-enhanced models have [Fe/H]=— 0.29. 
In the following discussion we will focus in detail 
on specific line indices which we will be using diag- 
nostically to determine iron content [Fe/H], total 
metallicity Z, and age. 

The left-hand panel of Figure [5] shows grids of 
H/3 vs. Fe5406 derived from our SSP integrated 
spectra for a range of metallicities and ages, as 
detailed in the figure caption. The figure includes 
both the scaled solar and a-enhanced grids and 
clearly demonstrates that the Fe5406 line is pre- 
dominantly a tracer of iron abundance, [Fe/H], 
since the 2 grids closely correspond along lines 
of very similar [Fe/H], and is generally insensitive 
to the degree of a enhancement. This is not the 
case for other Fe line indices we have tested, which 
also include other d ominant elements, as listed in 
Traeer et ail (|l998l their Table 2). Although our 
line indices are not converted onto the Lick system, 
this sensitivity of the Fe5406 line to Fe only is sim- 



3 See http://www.iac.es/galeria/vazdekis/index.html 
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ilar to that noted in other studies which do utilize 
i ndice s on the Lick sy s tem, e.g. 
(120051 ) and lKorn et all (l2005h . 



Lee fc Worthev 



The right-hand panel of Figure [5] shows similar 
grids for H/3 vs. [MgFe], where [MgFe] is defined as 
V< Fe > x Mg6 and <Fe>=±(Fe5270+Fe5335). 
It can clearly be seen that [MgFe] is sensitive to 
total metallicity, since the scaled solar and a- 
enhanced grids almost exactly correspond along 
lines of constant Z across a wide metallicity range. 
An important trend to note in this plot is the be- 
haviour of the H/3 line which, at fixed total metal- 
licity Z and age, is stronger in the a-enhanced 
models than the scaled solar ones. This behaviour 
of the H/3 lin e strength is qualitat ively similar to 
that found bv lTantalo et alj who used high 

resolution spectra to determine response functions 
which are used in conjunction with existing fit- 
ting functions to 'correct' solar scaled indices to 
the appropriate q-enhance ment (effectively, an up - 
date of the method used bv lTripicco fc Belllll995l ). 
However we note h ere that the recent work of 
Coelho et all (|2007l ). discussed in Section 12.1.11 



seems to display the opposite trend in the be- 
haviour of H/3 with respect to Z (see discussion 
below) . 

These 2 figures seen together (i.e. the 2 pan- 
els of Figure [6|) illustrate very well the difficulty 
in predicting ages from the H/3 line without also 
having other information, in particular, the degree 
of a enhancement. From our models it seems that 
the H/3-Fe5406 combination can be used to un- 
ambiguously determine age and [Fe/H], but yields 
no information about the degree of enhancement. 
Conversely, using H/3- [MgFe] alone could result in 
spurious ages unless the appropriate level of a en- 
hancement is determined by other means. Hence 
the diagnostic use of a single index-index diagram 
is not recommended. 

In Figure[7] (left-hand panel) we show the corre- 
sponding grids for H/3 vs. Mg6 which demonstrate 
the sensitivity of the Mg6 index to the degree of a 
enhancement (note that in our models this is fixed 
at [a/Fel=0. 4). Qualitatively similar behaviour is 
also seen bv ISchiavonl (|2007l ). However, this is a 
difficult grid to use diagnostically because the sen- 
sitivity to changing [Fe/H] and Z are impossible to 
disentangle from this grid alone and hence, in our 
later analysis, we will utilise the Fe5406-[MgFe] 
plane, illustrated in the right-hand panel of Fig- 



ure [7] In this plane, the lines of constant age are 
almost completely degenerate, however the scaled 
solar and a-enhanced grids are clearly separated. 
Using a combination of 3 grids - H/3 vs. Fe5406, 
H/3 vs. [MgFe], and Fe5406 vs. [MgFe] - it should 
be possible to disentangle age, total metallicity Z, 
[Fe/H] and the degree of a enhancement, as we 
demonstrate in Section l3~T1 

We compared the line strengths derived from 
our high reso l ution integrated spectra to those of 



Coelho et al.l (|2007[ ). a recent study which uses a 



method most similar to ours to produce integrated 
spectra for a range of SSPs (see Section 12.1.11 
for a comparison with their colors). The upper 
left panel of Figure [5] shows our scaled solar H/3 
vs. Fe5406 grid with the Coelho data superim- 
posed for [Fe/H]=Q.O (ages 3-12 Gyr). There is 
reasonably good agreement between their models 
and ours in the general locus of points along the 
[Fe/H]=0.0 line, however there are obvious offsets 
in H/3 which would cause a discrepancy in pre- 
dicted ages. In general, their scaled solar grid lies 
at higher H/3 values than ours and so for a fixed 
(observed) data point we would predict younger 
ages. The upper right hand panel of Figure [5] 
shows the full grid of H/3 vs. Fe5406 from the 
Coelho models (note, only 3 metallicities are cur- 
rently available) - this figure seems to show the 
opposite trend in the behaviour of Hp with re- 
spect to a-enhancement compared to our models, 
as mentioned above. 

We have no obvious explanation for this differ- 
ence in behaviour bu t it sho uld be noted that, al- 
though [CoeihoIetjiL ( 2007 ) use a method which 
is qualitatively similar to ours to produce inte- 
grated spectra, they are using different underly- 
ing stellar models and a different spectral library. 
In particular, we note here that the isochrones 
used in the Coelho models differ signficantly from 
the BaSTI ones when similar ages and metallici- 
ties are compared. In general, their RGBs are at 
cooler temperatures than the BaSTI ones, by up 
to ~ 200K. Also the temperature at the turn off 
(TO) displays significant differences in behaviour, 
e.g. at old ages, the BaSTI isochrones have hotter 
TOs, whilst at intermediate and young ages this 
trend is reversed, and the BaSTI TOs are cooler. 



Coelho et al.l ([20071 ) also include some low temper- 



ature spectra in their models, generated using the 
MARCS model atmospheres, in order to cover all 



S 



evolutionary stages - this enables them to com- 
pute broadband colors from their high resolution 
spectra, whilst we use our low resolution spectra 
for this. They also compute models for different 
values of total metallicity Z to ours, meaning that 
a comparison of grids (i.e. with a range of ages and 
metallicities) is very difficult - only their scaled 
solar, [Fe/H]=0.0 models are directly comparable 
with ours. 

For the interested reader, we also include a 
similar co mparison with the fitti ng-function based 

(|2003l ) - these are dis- 



models of iThomas et al 



played in the lower panels of Figure [8] By defini- 
tion, fitting-function based models are on the Lick 
system and so we have modified the tabulated H/3 
and Fe5406 values by 0.13 and 0.2 respectively, ac- 
cordin g to the prescription of iBruzual fc Charlotl 
(|2003l their Table 6), to make them directly com- 
parable with indices derived from flux calibrated 
spectra, as used in our work. 

We note here that all of our models presented 
so far in this paper have used the BaSTI non- 
overshooting models with a value for the Reimers 
mass loss parameter of r;=0.2. The effect of vary- 
ing r\ will be discussed in the next section. 

3. Comparisons with resolved stellar clus- 
ters 

Any stellar population synthesis model should 
be tested against resolved populations with inde- 
pendent estimates of chemical composition and 
age. Star clusters belonging to the Galaxy are 
the obvious choice, because their chemical compo- 
sition can be determined from spectroscopy, and 
theoretical isochrones matched to the observed 
color-magnitude diagram (CMD) provide an esti- 
mate of the age. Ideally, the analysis of integrated 
colors and spectra by means of population syn- 
thesis methods should provide a chemical compo- 
sition and age consistent with the in ferences from 



spectroscopy and the CMD (see, e.g.lGibson et al. 
19991: IVazdekis et al.ll200lUSchiavonf 2007. and ref- 
erences therein). 

To this purpose one can make use of the 
Schiavon et al.1 (|2005h library of integrated spec- 
tra of 40 Galactic globular clusters, plus the in- 
tegrat ed spectrum of the G alactic open cluster 
M67 (|Schiavon et al.l l2004ah . Each of the glob- 
ular cluster integrated spectra covers the range 



- 3350 - 6430A, with - 3.1 A (FWHM) resolu- 
tion, and one can potentially attempt fits to the 
whole sample of spectra, comparing the results de- 
rived from diagnostic lines to those derived from 
isochrone fitting to the observed CMD. There are 
however some serious issues that one has to take 
into account when interpreting the results of these 
comparisons. 

The first problem is the [Fe/H] estimates of 
the clusters. There are lar ge discrepanc i es be- 
twe en the widely employed Zinn fc West! (Il984h 
and Carretta fc Grattonl ( 1997 ) scales. For exam- 
ple, well studied objects like M3 and M5 show dif- 
ferences of ~0.3 dex between the two [Fe/Hl es- 
timates. A third [Fe/H] scale by iKraft fc Ivans! 



(|2003l ) based on Fell lines shows very large dif- 
ferences compared to the previous two sets of es- 
timates in the low metallicity regime. As an ex- 
ample, a well studied metal poor cluster like M68 
displays a ~ 0.4 dex range of [Fe/H] values when 
these three different [Fe/H] scales are employed. 
Taking a mean value of the various determinations 
for each cluster does not make much sense, given 
that the differences among the authors are mainly 
systematic. In addition, [a/Fe] spectroscopic esti- 
mates do not exist for many clusters. 

The second issue is the presence of the well 
known CN and ONa anticorrelations in the metal 
mixture of st ars within a single cl uster (see, e.g., 



the review by iGratton et alj 120041 . and references 
therein). There is now convincing evidence that 
C, N, O, Na - and sometimes also Mg and Al - 
display a pattern of abundance variations super- 
imposed onto a normal a-enhanced heavy element 
distribution. Negative variations of C and O are 
accompanied by increased N and Na abundances. 
There is also now general agreement that these 
abundance patterns are of primordial origin. This 
means that the integrated spectrum of a cluster 
is composed of individual spectra with a range of 
values for the very important CNONa elements, 
whereas the Fe abundance is essentially constant 
for all stars. It is therefore very difficult to in- 
terpret the strength of all line indices affected by 
these 4 elements, given that in a cluster the exact 
proportions of stars with a certain set of CNONa 
abundances, and their location along the CMD, 
are known only for small samples of objects. 

The third issue is related to the morphology 
of the horizontal branch (HB). As is well known 
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(see, e.g. Schiavon et al. l2004bh . the presence of 
blue HB stars affects the Balmer lines, which can 
produce spuriously young ages for clusters with 
an extended blue HB if this is not included in 
the theoretical modelling. A detailed synthetic 
modelling of the HB color extension for each of 
the observed clusters is therefore in principle nec- 
essary. Whilst this is possible in principle, it is 
difficult to implement in practice and would be 
computationally very time consuming to model all 
the different possibilities. This approach also as- 
sumes that one knows, a priori, the appropriate 
morphology required - this of course is not possi- 
ble for unresolved stellar clusters. The color of 
the HB is determined by the mass loss history 
of RGB stars when the initial chemical compo- 
sition is fixed. In the theoretical isochrones, this 
is modelled using the Reimers mass loss law, and 
the mass loss parameter 77 is a free parameter in 
the models, which is set at some fixed value for 
each isochrone. This means that for any particular 
isochrone, all the stars on the HB have essentially 
the same mass and hence start their evolution from 
the Zero Age HB all with the same color. For any 
fixed r\ value, higher metallicity SSPs have a red 
clump of stars along the HB. The distribution of 
objects along the HB becomes progressively bluer 
when the metallicity decreases, provided that 77 is 
unchanged. In Galactic globular clusters however, 
the mean color and color extension of the HB is 
often not correlated to the cluster metallicity, and 
additional chemical or environmental parameters 
that directly or indirectly affect the mass distri- 
bution along the HB, must be at play (see e.g . 
Caloi fc D'Antonall2005l : lRecio-Blanco et al.ll2006L 



library by ISchiavon et al.l (|2005l ) to our high res- 
olution spectra calculated with both i] — 0.2 and 
j] = 0.4. We have measured indices on the GC 
spectra in exactly the same way a s for our model 
spectr a. We note here that the ISchiavon et al 



and references therein). If 77 is fixed at a higher 
value in the models, the blue extension of the HB 
becomes more extreme at low metallicity, partic- 
ularly for old ages (see discussion below). 

A final issue is the existence of blue stragglers 
in the CMDs of Galactic globular clusters, i.e. a 
plume of stars brighter and generally bluer than 
the TO, thought to be created through mass trans- 
fer in a binary system or through direct collision 
of 2 stars. Integrated spectra that include only 
the contribution of single star evolution miss this 
component, causing potential biases in their age 
estimates. 

Keeping in mind all these limitations, we first 
compare the whole globular cluster (GC) spectral 



([20051 ) data set includes multiple observations for 
several clusters - in the following plots all the data 
points are shown. To minimize the effect of the 
CNONa anticorrelations on the comparison, we 
consider first the Fe5406— [MgFe] plane to esti- 
mate [Fe/H] (which is constant among all stars 
in a cluster) and the degree of a-enhancement 
from our theoretical spectra. As seen in sec- 
tion [272H1 [MgFe] is sensitive to the global metal- 
licity of the stellar population, but not to the ex- 
act distribution of the metals, and therefore the 
Fe5406— [MgFe] diagram - insensitive to age for 
old stellar populations - depends on the relation- 
ship between the global metallicity Z and the iron 
content in the population. Given that at fixed 
[Fe/H] an a-enhanced mixture has a larger Z com- 
pared to a scaled-solar one, the position of a clus- 
ter in this diagram is related to the degree of a- 
enhancement in the spectra of its stars. A possible 
complication is that the CNONa anticorrelation 
modifies the baseline a-enhanced metal mixture 
in a fraction of stars within the cluster. However, 
if the sum of C+N+O (which makes up a sub- 
stantial fraction of the total metals) is unaffected 
by the abundance anomalies (as seems to be the 
case within current observational errors, see, e.g. 
ICarretta et al. 2005 ) the relation between [Fe/H] 
and Z is unchanged, and the Fe5406— [MgFe] di- 
agram still provides an estimate of the degree of 
a-enhancement in the baseline chemical composi- 
tion of the cluster as a whole. 

Figure [9] shows all the Schiavon data points in 
the Fe5406-[MgFe] plane, with lines of 14 Gyr 
(constant) age derived from our scaled solar and a- 
enhanced spectra, in the metallicity range — 1.8 < 
[Fe/H] < +0.06. In general, the data points fall 
exactly between the scaled solar and a-enhanced 
lines indicating that most of them are probably a- 
enhanced to some extent, as expected. The group 
of points at the top of the grid represents 2 clusters 
only - NGC6528 (for which there are 6 observa- 
tions) and NGC 6553, which are the 2 most metal 
rich clusters in the data set at [Fe/H]^ —0.2. The 
4 points at the bottom left of the grid represent 
2 observations each of NGC 2298 and NGC 7078, 
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the 2 most metal poor clusters at [Fe/H]~ —2.0. 
Hence we maintain that this grid is a good first 
indicator of [Fe/H], and also provides an estimate 
of the degree of a enhancement. 

Once [Fe/H] and the degree of a-enhancement 
are established we can then use the H/3-Fe5406 
plane to estimate the cluster ages - Figure [10] 
shows the same GC data points on our a-enhanced 
grid, calculated using 77=0. 2. Whilst the interme- 
diate metallicity clusters generally span the 10- 
14 Gyr age range (within the errors), it can be 
seen that the more metal poor clusters all seem to 
scatter to higher ages. This is most likely to be 
due to the blue extent of their HBs, as discussed 
above. To illustrate the potential effect of having 
extended blue HBs due to increased mass loss on 
the RGB, we have overplotted the 14 Gyr line cal- 
culated using 77=0. 4 - this clearly demonstrates 
the large increase in the strength of the H/3 line 
which can lead to spuriously young ages. We note 
that the trend of this 77=0. 4 line, i.e. the rise to a 
peak in H/3 at some intermediate metallicity and 
then a decrease towards lower metallicities, is not 
unexpected. For the most metal poor old clusters, 
the horizontal branch can develop a very long blue 
tail which extends down to magnitudes almost as 
faint as the turn off in optical bands. Thus al- 
though the individual stars still exhibit strong H/3 
due to their very hot temperatures, their contri- 
bution to the overall continuum flux decreases at 
the wavelength of H/3 (which is in the middle of 
the -B-band) at the lowest metallicities. 

3.1. Case study: 47 Tuc 

We have performed a more detailed com- 
parison for the case of 47 Tuc, which is one 



the ISchiavon et"aH (|2005h 



of the clusters 
database. This c luster has a reasonably we ll es- 
tablished [Fe/H]. ICarretta fc Grattonl Jjjjjl find 
[Fe/H]=-0.70 ± 0.07. IZinn fc Westl dl984h give 
[Fe/H]=-0.71 ± 0.08, while iKraft fc Ivand (|2003h 
obtain between —0.78 and —0.88 depending on 
the assumed T e ff scale, based on the equivalent 
widths of Fell lines measured from high-resolution 
spectra of gian ts. More recent hig h resolution 
spectroscopy bv lCarretta et aL (2004) gives an es- 
timate consist ent with that of Carretta fc Gratton 
<ll997h . while iKoch fc McWilliaml (|2008h obtain 
[Fe/H]=-0.76 ± 0.01 ± 0.04. Measurements of a- 



iKoch fc McWilliaml (|2008h give a enhancement, 
[a/Fe], of the order of ~ 0.3 — 0.4 dex on average. 

Figur e [TT1 displays a fit to the BV data by 
Stetson! ((2000) with the BaSTI a-enhanced, 77=0.2 
isochrone and ZAHB, with [Fe/H]=— 0.7 and an 
age of 12 Gyr. A more formal determination of the 
age and associated errors has bee n obtained using 
the Ay met hod described in, e.g. JSalaris fc Weissl 
(| 19971 . 119981 ) . This method employs the difference 
in V magnitude between the Zero Age HB (ZAHB) 
and the turn off (TO), designated Ay, as an age 
indicator. This quantity turns out to be weakly 
sensitive to the cluster metallicity, so that errors 
on the age estimate are dominated by errors on the 
determination of the ZAHB and TO magn itudes 
(see the discussion in ISalaris fc Weissl Il997l) . We 
obtain for 47 Tuc an age t = 12±1 Gyr where in 
the error budget we have included the small con- 
tribution of a ±0.2 dex uncertainty around the 
reference [Fe/H] =—0.7. 

A comparison of the 47 Tuc line indices with 
those derived from the theoretical spectra (com- 
puted with 77=0.2) in the Fe5406-[MgFe] and 
H/3— Fe5406 diagrams is included in Figures [9] and 
fTUl We have estimated representative error bars 
for 47 Tuc, which will be approximately the same 
for all the GC data points. These include con- 
tributions from the measurement of the line in- 
dices on the observational data as given by the 
LECTOR program, an estimate of the intrinsic 
observational error (there is only one observation 
of 47 Tuc), and a small contribution due to the 
difference in resolution between the data and the 
models. Within the error bars of the observed in- 
dices, our BaSTI integrated spectra yield [Fe/H], 
[a/Fe] and age estimates in agreement with the 
results from high resolution spectroscopy of indi- 
vidual stars and from the CMD age. 

We have then checked the extent of possible 
biases introduced by the theoretical representa- 
tion of the cluster HB, which differs from the 
observed one. To this purpose we have consid- 
ered the Fe5406, [MgFe] and H/3 values obtained 
from a 14 Gyr old, 77=0.2, [Fe/H]=-0.7 theoret- 
ical spectrum. We have then computed an inte- 
grated spectrum for the same age and [Fe/H], but 
co nsidering a synthet ic HB (obtained as described 



element abundances bv lCarretta et alj (|2004[ ) and 



in IS alar is et al.l 120071 ) that precisely matches the 
color extension observed in 47 Tuc. The result- 
ing Fe5406, [MgFe] and H/3 values are practically 
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unchanged compared to the 77=0.2 case - the in- 
crease in H/3 is < 0.04, corresponding to an age 
difference of < 0.5 Gyr, and well within the ob- 
servational error, whilst Fe5406 and [MgFe] are 
negligibly changed. 

As a second check we have estimated the ef- 
fect of the blue stragglers present i n the cluster 



sis. The contribution of blue stragglers is excluded 
from the spectrum we employ in the comparisons 
that follow. 

For this cluster the s pectr os copic invest i gation s 



CMP . According to the estimates bv lPiotto et al 
([20021 ) . the ratio of blue stragglers to the num- 
ber of HB stars is about 0.15, one of the lowest 
fractions within the Galactic globular cluster sys- 
tem. To determine an upper limit of their effect 
on the indices used in our analysis, we have added 
to the 14 Gyr old [Fe/H]=-0.7, 77=0.2 integrated 
spectrum, the contribution of a clump of Zero Age 
MS stars 1.5 mag brighter than the isochrone TO, 
which is the approximate position of the bluest 
and brightest blu e stragglers in the BV CMD by 
Sills et all (l2000h . Their number is assumed to be 
15% of the number of HB stars predicted by the 
models. All line indices measured on the result- 
ing integrated spectrum are practically unchanged 
compared to the reference values, with differences 
of no more than 0.01 in any of the diagnostic in- 
dices used here. 

3.2. Case study: M67 

As a second detailed t est, we have co n sidered 



the open cluster M67. ISchiavon et al.1 (|2004al ) 
have created a representative integrated spectrum 
for M67, built up from observations of spectra of 
individual stars in the cluster (about 90 objects) 
that sample as well as possible the CMD from 
the lower MS to the tip of the RGB and the He- 
burning (red clump) phase. Some local field stars 
are also included to better sample some regions 
of the cluster C MD not well represente d by the 
M67 data alone. Schiavon et al.1 ( 2004ah assigned 
a mass to each of the sample stars by using inputs 
from theoretical isochrones, and then the individ- 
ual spectra were coadded with weights given by an 
adopted IMF (they employ a Salpeter IMF). The 
contribution of objects with masses below ^0.85 
M© (below the f aintest objects whose spectrum 
is included in the ISchiavon et al. 2004a! analysis) 
has been included by Schiavonl^OO? ) using in- 
puts from theoretical models. The corrections that 
need to be applied to line indices measured on 
their published M67 spectrum to account for this 
(their Table 11) have been included in our analy- 



i<or tins cluster tne s pectr os copic inv estigation; 
bv iTautvaisiene et al.l (|2000h . lYong et all (120051) 



Randich et al.l (|2006T) and the analysis by Taylor 
((2007j), all give [Fe/H]~ 0.0. The best age esti- 
mate obtained by fitting our [Fe/Hl=+0.06 BaSTI 



non-overshooting isochrones to ISandauist ( 2004 ) 
photometry is £=4.0 Gyr. The inclusion of over- 
shooting has only a marginal effect at this age, 
increasing t to 4.2 Gyr. The quality of the fit is 
very similar, with a small preference for the non- 
overshooting case, which we display in Figure 1121 
A comparison of the M67 line indices with those 
derived from our scaled solar non-overshooting 
grid (computed with 77=0.2, but at the age of M67 
the choice of 77 is irrelevant) in the H/3— <Fe> di- 
agram is displayed in Figures fT3l We used <Fe> 
instead of the Fe5406 line here because the wave- 
length coverage of the M67 spectrum is slightly 
less than for the GC sample, cutting off at 5396A. 
We have not assigned error bars to the data point 
on this plot as it not trivial to quantify them - 
they must however be at least as large as the error 
bars estimated for the GC data, and some other 
potential contributory factors are discussed below. 
Whilst the age inferred from this diagram is con- 
sistent with the CMD age, the predicted [Fe/H] 
is lower than the spectroscopic value by about 
0.2 dex. There are some important factors which 
we think may contribute to this discrepancy which 
we address in the following discussion. 

One important point to consider is that M67 
spectrum does not come from observations of the 
cluster integrated light (as is the case for the GC 
data set), rather it has been built in the same 
way as theoretical integrated spectra are built, 
starting from an isochrone. Each point along 
the isochrone has been assigned an appropriate 
spectrum, and the contributions of the individual 
spectra weighted by the chosen IMF. The BaSTI 
isochrones include more than 2100 points to sam- 
ple masses above 0.85M Q - a very fine coverage 
- whereas the M67 spectrum has been computed 
using a much smaller number of individual rep- 
resentative stars. To get a handle on the effect 
introduced by the sampling used for the observed 
spectrum, we have computed two integrated spec- 
tra for the same age of 4.0 Gyr ([Fe/H]=0.06, 
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scaled solar non-overshooting models). The first 
one has been obtained as usual, whereas for the 
second spectrum the resolution of the isochrone 
for masses above 0.85M Q has been degraded to in- 
clude only 90 points. The spacing of these points 
has been chosen to reproduce as well as possible 
the V magnitu d e distri bution of the stars used by 
Schiavon et al.1 ( 2004ah (we have assumed a dis- 
tance modulus (m — M)o=9.67, as obtained from 
the isochrone fitting) . We find a non-negligible dif- 
ference between line indices computed with these 
two different isochrone samplings. The value of 
the difference between the fine sampling and the 
coarse sampling spectrum has then to be added 
to the measured values for M67, in order to have 
a more consistent comparison with the theoretical 
grid. The direction in which the data point moves 
in the the H/3— <Fe> plane after applying this 
correction is indicated on Figure [13] by the arrow. 
The H/3 value decreases, implying a slightly older 
age, although this may not be significant within 
the observational errors. Interestingly, all the Fe 
indicators move to significantly higher values, in- 
creasing the inferred [Fe/H] by as much as 0.2 dex 
and putting the predicted [Fe/H] more in line with 
the spectroscopic value. 

This test highlights some subtle issues affecting 
the creation of SSPs which can have non-trivial 
effects on their derived integrated properties and 
spectra. We remind the reader that the method 
of populating an isochrone to create an SSP used 
here (and in other population synthesis studies) 
is an analytical one, in which the isochrone is ef- 
fectively treated as a continuous function. This 
method is strictly correct only in the limit of an 
infinite number of objects in the stellar system be- 
ing modelled, but is the appropriate method to 
use when modelling an unresolved population. In 
practice, an isochrone consists of a discrete num- 
ber of evolutionary points (EPs), all of which are 
assigned a spectrum when the analytical method 
is used. However, the number of defined EPs, re- 
sulting in coarse or fine sampling of the isochrone, 
can make a difference to the resultant integrated 
properties, as demonstrated above. 

A further refinement in the modelling of the in- 
tegrated properties of observed stellar clusters can 
be made by populating the isochrone with a dis- 
crete number of stars appropriate to the type of 
cluster (using an IMF, as before). The Web tool 



provided on the BaSTI database which creates 
synthetic CMDs for resolved populations (SYN- 
THETIC MAN) populates isochrones in this way, 
with the number of stars in the simulation being 
determined by the user. When the number of stars 
in a particular simulation is not large enough to 
populate every EP on the isochrone (which is al- 
ready true for typical 10 5 — 10 6 M Q Galactic glob- 
ular clusters) statistical fluctuations of star counts 
will arise. This means that an ensemble of SSPs 
all with the same total mass, age and initial chem- 
ical composition will display a range of integrated 
magnitudes and colors, that increases (for a fixed 
value of Mt) when moving towards longer wave- 
lengths, which are dominated by the shorter-lived 
(hence more sparsely populated) RGB and AGB 
phases. Preliminary tests by us have shown that 
there are also non-negligible effects on integrated 
spectra produced in this way and that a signif- 
icant scatter is seen in derived line indices when 
the number of stars in the simulations corresponds 
to a total mass less than ~ 2 x 10 5 M Q (i.e. roughly 
the size of a typical GC and larger than most 
open clusters). In particular, one significant fac- 
tor seems to be the exact location of the brightest 
star populating the RGB, which can easily vary 
by several tenths of a magnitude when modelling 
open clusters like M67. 

All these issues, and particularly the effects of 
statistical fluctuations on integrated spectra and 
derived line indices, will be explored further in a 
subsequent paper. 

4. Summary 

In this 4th paper in the BaSTI series we have 
presented a new set of low and high resolution syn- 
thetic spectra based on the BaSTI stellar models, 
covering a large range of SSPs, for both scaled so- 
lar and a enhanced metal mixtures. This enables 
a completely consistent study of the photometric 
and spectroscopic properties of both resolved and 
unresolved stellar populations, and allows us to 
make detailed tests on specific factors which can 
affect their integrated properties. 

• We have produced low resolution spectra 
suitable for deriving broadband magnitudes 
and colors in any photometric system and 
ensured that the derived colors are fully con- 
sistent with those predicted from the BaSTI 
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isochrones. 

• We have produced high resolution spectra 
suitable for studying the behaviour of line in- 
dices and tested them against a large sample 
of Galactic globular clusters, making partic- 
ular note of the effect of horizontal branch 
morphology. 

• We have tested the global consistency of the 
BaSTI models by making detailed compar- 
isons between ages and metallicities derived 
from isochrone fitting and from line index 
strengths for the Galactic globular cluster 
47 Tuc and the open cluster M67. 

• We find non-negligible effects on derived line 
indices caused by statistical fluctuations, 
which are a result of the specific method 
used to populate an isochrone and assign 
appropriate spectra to individual stars. 
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A. Appendix 



Upon publication of this paper, the official BaSTI website (http : //193 . 204 . 1 . 62/index 7html[ ) will be 



expanded to include integrated spectra (low and high resolution) , magnitudes and mass-to-light (ML) ratios 
for SSPs (normalized to IMq of stars formed during the burst of star formation) covering the full range 
of ages, chemical compositions, choices of mass loss and mix ing included in the d atabase. A full summary 
of the parameter space covered by BaSTI can be found in ICordier et al. Integrated magnitudes 



and ML ratios will be provided for UBVRIJHKL filters. Additional photometric systems can be easily 
accounted for by convolving the BaSTI low resolution integrated spectra with the appropriate filter response 
curves, and using an appropriate normalizing spectrum (a Vega synthetic spectrum is also included in 
BaSTI). ML-ratio files also tabulate the amount of mass locked into the various types of remnants for each 
SSP, to help determine ML-ratios in any arbitrary photometric system. We remind the reader that the 
ML-ratio of a stellar population is customarily defined as the ratio between the integrated luminosity in a 
given passband (in solar units) divided by the integrated stellar mass (in solar units) contained in stars, 
including all remnants like white dwarfs (WDs), neutron stars (NSs) and black holes (BHs). To compute 
the solar luminosity in the UBVRIJHKL we have adopted the following set of solar absolute magnitudes: 
M v = 5.66, M B = 5.49, M v = 4.82, Mr = 4.45, Mj = 4.10, Mj = 3.68, M H = 3.31, M K = 3.28, M L = 3.26. 
The mass contained in WD, NS and BH remnants has been computed by assuming an initial-final-mass 
(IFMR) relation in addition to the IMF. For WDs with a carbon-oxygen core we employed the IFMR 
obtained from our AGB modelling. The upper initial mass limit for the production of carbon-oxygen WDs 
(the so-called M up ) is also obtained from our stellar models. For initial masses between M up and 10 M Q we 
considered typical oxygen-neon core WDs with masses of 1.3M Q . Initial masses between 10 and 25 M Q are 
assumed to give rise to 1.4M Q NSs, whereas for progenitors above 25M Q we consider a final state as black 
holes with mass equal to 1/3 of the initial mass. 
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Fig. 1. — Sample low resolution (upper panel) and 
high resolution (lower panel) integrated spectra 
from scaled solar, Z=0.0198 (the solar value) mod- 
els. Ages displayed are 3 and 10 Gyr, which are 
the upper and lower spectra respectively on each 
plot. 



Fig. 3. — Compari son of BaSTI c o lors (s olid lines) 
with those from Coelho et al ] (|2007h (dashed 
lines). The BaSTI colors are from our scaled so- 
lar [Fe/H]=+0.06 and a-enhanced [Fe/H]=+0.05 
models. The Coelho models are both for 
[Fe/H]=0.0. 



l 1 f— 




Fig. 2.— SSP colors from BaSTI models. Solid 
lines are the scaled solar models for [Fc/H]=+0.06, 
—0.66 and —1.27 (upper, middle and lower lines 
respectively on each plot). Dashed lines are the 
a-enhanced models for [Fe/H]=+0.05, —0.70 and 
-1.31. 
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Fig. 4. — Comparison of MARCS (solid lines) 
and Munari (dashed lines) models (scaled solar, 
[Fe/H]=0.0) for H/3, Mg6 and Fe5406, at effective 
temperatures T e// =4000K (red), 4500K (blue), 
5000K (cyan). 
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Fig. 5. — Percentage of flux in stars with T e ft < 
3500K for 2 representative SSPs. Bandpasses for 
H/3, Mgb and Fe5406 are marked with the dotted 
lines. 




Fig. 7.— H/3 vs. Mgb (left panel) and Fe5406 vs. 
[MgFe] (right panel). Models and symbols are as 
for Figure [6] 
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Fig. 6.— Grids of H/3 vs. Fe5406 (left panel) and 
[MgFe] (right panel) measured on our high reso- 
lution spectra as described in the text. The solid 
line/circles are the scaled solar grid, the dashed 
line/open squares are the a-enhanced grid. Scaled 
solar models are for 2"=0.0003, 0.0006, 0.001, 
0.002, 0.004, 0.008, 0.0198, 0.04 (corresponding 
to [Fe/H]=-1.79, -1.49, -1.27, -0.96, -0.66, 
—0.35, +0.06, +0.40) and a-enhanced models are 
for Z=0.0006, 0.002, 0.004, 0.008, 0.0198, 0.04 
(corresponding to [Fe/H]=-1.84, -1.31, -1.01, 
-0.70, -0.29, +0.05). Metallicity increases from 
left to right. Ages (increasing from top to bottom) 
are 1.25, 3, 6, 8, 10 and 14 Gyr. 



Fig. 8. — Compari son of our H / 3 vs. Fe5406 
grids with those of ICoelho et all (|2007t ) (upper 
panels) and iThomas et all (|2003f ) (lower panels). 
The left panels show the Coelho (upper) and 
Thomas (lower) [Fe/H]=0.0, scaled solar models 
for ages 3-12 Gyr superimposed on our scaled 
solar, H/3 vs. Fe5406 grid. Upper right shows 
the Coelho grids for scaled solar (solid) and a- 
enhanced (dashed lines) for [Fe/H]= —0.5, 0.0 
and +0.2. Lower right shows TMB03 scaled 
solar grids for [Fe/H]=-0.33, 0.0, 0.35 (solid) 
and a-enhanced ([a/Fe]=0.3), [Z/H]=0.0, 0.35, 
0.67, corresponding to [Fe/H]=-0.28, 0.06, 0.38 
(dashed). 
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Fig. 9. — Fe5406 vs. [MgFe] from our scaled solar 
and a-enhanced 14 Gyr spectra (solid and dashed 
lines, respectively) joined at approximately equal 
[Fe/H] (dotted lines) for [Fe/H] - +0.06, -0.3, 
—0.7, —1.0, —1.3, —1.8 (decreasing from top to 
bottom). Open circles are the Galact ic globular 



(2005|)-47Tuc is 



cluster data from lSchiavon et al. 
represented by the solid square, with an estimated 
typical error bar. 



Fig. 10. — H/3 vs. Fe5406 grid for our a- 
enhanced models ([Fe/H]=-1.84, -1.31, -1.01, 
-0.70, -0.29 +0.0 5) with GGC data from 
Schiavon et al.l ( 20051 ) overplotted - symbols as 



for Figure [U The grid shown is calculated using 
mass loss parameter rj—0.2. The dashed line is for 
the 14 Gyr models only and was calculated using 
7^=0.4, this illustrates the effect of a very extended 
blue horizontal branch on the predicted H/3. 
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Fig. 1 1 . — Is ochrone fit to 47 Tuc (data from 
IStetsonl l2000h . The best fit to a 12 Gyr, 
[Fe/H]=— 0.7, a-enhanced isochrone yields a dis- 
tance modulus of (to — Af)o=13.28, for an as- 
sumed reddening of E (B-V) =0.02, in agreement 
with the empi r ical m ain-sequence fitting result of 
Percival et all (|2002t) . 
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Fig. 12.— Is ochrone fit to M67 (data from 
Sandquist 2004 ). The best fit to a scaled solar, 
non-overshooting, [Fe/H]=0.06 isochrone, yields 
an age of 4.0 Gyr and a distance modulus of 
(to — M)o=9.67, for an assumed reddening of 
i?(B-V)=0.02, in agreement with the empirical 
main- sequence fitting result of iPercival &: Salarid 
(|2003l) . 
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Fig. 13. — Scaled solar grid for H/3 vs. <Fe> with 
the M67 data point, calculated as described in the 
text. The arrow shows the approximate direction 
and distance the point would move if the effects of 
statistical fluctuations were taken into account. 
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